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Fragmentation functions (FFs) describe the formation of final state particles from a par- 
tonic initial state. Precise knowledge of these functions is a key ingredient in accessing 
quantities such as the nucleon spin structure in semi-inclusive deep-inelastic scattering 
and proton proton collisions. However, fragmentation functions can currently not be de- 
termined from first principles Quantum Chromodynamics and have to be extracted from 
experimental data. The Belle experiment at KEK, Japan, provides a large data sample 
for high precision measurements on e^e~ annihilations allowing for first-time or more pre- 
cise extractions of fragmentation functions. Analyses for extractions of spin-independent 
(unpolarized FFs) as well as spin-dependent fragmentation functions (interference FFs) at 
Belle are presented. 



1 Precision measurement of charged pion and kaon mul- 
tiplicities 

1.1 Motivation and outline 

A number of polarized Semi-Inclusive Deep-Inelastic Scattering (SIDIS) and polarized proton- 
proton scattering (pp) experiments are measuring quark and gluon distributions via QCD anal- 
ysis. Unpolarized hadron fragmentation functions are input quantities for such analyses and 
currently limit the accuracy with which e.g. parton helicity distributions can be determined. 
Two recent studies have extracted hadron fragmentation functions from datasets of hadron pro- 
duction in e+e~ [T, and from a combination of e+e~, SIDIS and pp datasets [T. In both studies, 
the authors emphasize the need for a precision hadron multiplicity measurement from e+e~ an- 
nihilation data at low center-of-mass energies, compared to the bulk of existing datasets taken 
at the Large Electron-Positron Collider (LEP) at CERN. This precision dataset is expected to 
improve in particular the knowledge of the gluon fragmentation functions. 

Presented in this section is a precision measurement of hadron multiplicities on about 220 x 
10^ e+e~ annihilation events taken with the Belle experiment at KEK, Japan, at a center- 
of-mass energy of 10.52 GeV, 60 MeV below the T(46') resonance. The multiplicities are 
measured in dependence of 2: = 2Ehad/y/s = Ehad/ Eparton-, the hadron energy relative to half 
of the center-of-mass energy, on an interval 0.2 < z < 1.0. The measured hadrons are produced 
in reactions e+e~ qq^ where q = {u^ 5, c}. 
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a) b) 

Figure 1: Fits of experimental invariant mass distributions rriD* — m]jo from decays D* 
TT -\- ^ TT + (Ktt). Figures a) and b) show distributions containing negatively-charged kaon 
candidate tracks with laboratory frame momentum 1.4 < piab < 1-6 GeV/c and laboratory 
frame azimuthal angle 77.9 < ^lab < 89.0 degrees. Figure b) shows all contributions to Figure 
a) where the kaon candidate additionally passes PID likelihood cuts to select pions. The PID 
misidentification of kaons as pions can be extracted from the ratio of the hatched areas, yielding 
the probability p^K-^n-) = 0.111 ± 0.004. 

1.2 Experimental-data-based calibration of the Belle particle identi- 
fication 

The Belle detector contains several subsystems which allow for particle identification (PID) 
by imposing cuts on likelihood values extracted from measurements of these subsystems. The 
likelihood cut selections yield fairly accurate PID but need to be calibrated in the context of 
a high precision measurement, such that measured hadron yields can be corrected for particle 
misidentification. Particle misidentification changes hadron yields up to 10% for pions and up 
to 20% for kaons, depending on hadron momentum. The correction is performed through an 
unfolding technique based on inverse 5x5 PID probability matrices. 

PID probability matrices are obtained from analyzing decays of particles in which the species 
i of the decay products can be determined through purely kinematic means. In such a sample, 
additional cuts on track PID likelihood variables selecting a species j define a subsample of 
tracks. Comparing the number of cut-selected tracks to the total number of tracks in the sample 
allows to extract PID probabilities p(^i^jy Complete 5x5 matrices of PID probabilities are 
extracted from experimental data for species {i, j} = {e^ jj^^tt^ K^p} by analyzing decays of D* ^ 
A and J/i/j particles. This data-driven PID calibration avoids the dependence on the modeling 
of Belle PID detectors in GEANT [3 . As an example. Figure [T] displays fits of invariant mass 
distributions tud* — m]jo from decays ^ tt + ^ tt + (Kir). The kaon candidates in 
contributions to Figure b) fulfill an additional PID likelihood cut to select pions. The ratio of 
the hatched peak areas is equal to the PID probability P^k-^tt-)- 

For kinematic regions not accessible to this calibration method, an extrapolation algorithm 
is used to obtain PID calibration information. The PID probability matrices from all kinematic 
areas are inverted and then applied to the measured yields 

meas — 
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Figure 2: a) Measurement of hadron multiplicities: Preliminary negatively-charged pion and 
kaon multiplicities including statistical and systematic uncertainties. An additional 1.4% nor- 
malization uncertainty is not shown, b) Measurement of the interference fragmentation func- 
tion: Azimuthal modulations al2 of two-pion yields for a 9 x 9 zi, ^2 binning, as a function of 
Z2 for all z\ bins. Indices 1 and 2 refer to the respective reaction hemisphere. 



obtain PID-corrected pion and kaon yields NpiD-corr- All uncertainties of the extracted PID 
probabilities are propagated through the inversion process using a Monte Carlo technique and 
are assigned to the PID-corrected yields as systematic uncertainties. 



1.3 Other systematic corrections, preliminary results and outlook 

In addition to particle misidentification, the measured experimental data yields are also cor- 
rected for sample purity, kinematical smearing, dec ay- in- flight, detector interaction/shower par- 
ticles, detector and tracking efficiencies and analysis acceptance. The presence of initial state 
radiation is accounted for as well. Figure [2] a) shows preliminary negatively charged pion and 
kaon multiplicities including statistical and systematic uncertainties from about 220 x 10^ e+e~ 
annihilations. In Figure [3] a) and b) the relative size of preliminary statistical and systematic 
uncertainties is given. Final systematic uncertainties will likely remain below 2% (3%) for tt 
(K) spectra for z < 0.6, and will increase with z up to 14% (50%) for tt (K) spectra through 
z ~ 0.9, respectively. Currently, the last systematic studies are being finished and timely 
journal submission for publication is expected. 



2 Measurement of the interference fragmentation func- 
tion for charged pion pairs at Belle 

Spin-dependent, chiral-odd fragmentation functions can be used to extract transverse spin quark 
distributions (so called transversity distributions) in the nucleon from polarized SIDIS and 
pp scattering experiments. One candidate for such a chiral-odd function is the interference 
fragmentation function (IFF). First proposed by Collins, Heppelman and Ladinsky [J , the 
IFF describes the fragmentation of a polarized quark into two hadrons correlated via partial 
wave interference. The corresponding observable is an azimuthal dependence in the production 
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a) b) 

Figure 3: Relative uncertainties of preliminary negatively-charged pion (a) and kaon (b) mul- 
tiplicities. An additional 1.4% normalization uncertainty is not shown. 

cross-section of hadron pairs. At Belle, the product of two IFFs can be measured by identifying 
one hadron pair in either hemisphere in a two-jet event e+e~ qq ^ {hh){hh)X . The IFF 
does not depend on transverse momenta and therefore its factorization and evolution can be 
described in a collinear approach. This makes the IFF an attractive alternative compared to 
the usually used Collins fragmentation function (previously measured for pions at Belle, cf. 
References and [6]) for the extraction of parton transversity distributions. 

Azimuthal correlations between the planes spanned by two hadron pairs and the interac- 
tion plane defined by the incoming lepton momenta and the quark- ant iquark momentum axis 
(approximated by the thrust axis) are extracted. The raw azimuthal yields are normalized to 
the average hadron pair yield and fitted. This approach is not susceptible to QCD radiative 
effects and the cos{(j)i + 02 ) azimuthal modulation ai2 is directly proportional to the product 
of two IFF Hf{z,m), where z = E pair / E parton represents the fractional energy of the hadron 
pair and m its invariant mass. Analysis results on a dataset of 672 fh~^ containing 711 x 10^ 
7r+7r~ pairs collected at the T(4S) resonance have been published in Reference 0. Figure [2] 
shows azimuthal modulations rising with the fractional energy of the hadron pair to significant 
values of up to 10% for highest z. 

The Belle Collaboration is currently also pursuing measurements of the kaon Collins FF, 
the kaon and tt^ IFFs and the di-hadron FF. In addition, measurements are being prepared to 
extract the kt dependence of the unpolarized and the di-hadron FFs. 
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